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ABSTRACT 


A series of measurements have been carried out with the purpose 
of quantitatively determining the characteristics of a commercial 
6.8 cm access superconducting rock magnetometer located in the magnetic 
properties laboratory at the Goddard Space Flight Center. The measure- 
ments have shown that although a considerable improvement in measurement 
speed and signal to noise ratios can be obtained using such an instrument, 
a number of precautions are necessary to obtain accuracies comparable 
with more conventional magnetometers . These include careful calibration 
of the sensor outputs, optimum positioning of the sample within the 
detection region and quantitatively establishing the degree of cross- 
coupling between the detector coils. In order to examine 
the uniformity of response for each detector, the responses have been 
mapped as a function of position, using a small dipole. The response 
variations were found to be less than 5% within an optimally positioned 
cylinder 2.5 cm in diameter and 4 cm in length. With the application of 
suitable correction terms, the magnetization of a 2.5 cm diameter core 
sample can be determined to an angular accuracy of better than half a 
degree, using a single insertion of the sample. For smaller diameter 
access magnetometers, spatial variation of the coil response characteristics 
will be more significant requiring repeated measurements of inhomogeneous 
samples in different orientations to average out the resulting errors. 



INTRODUCTION 


Superconducting rock magnetometers have, in the last few years, 
come from being something of a scientific curiosity, to a standard 
instrument used in paleomagnetic and rock magnetic research. 

The reason for this is their capability of measuring weakly 
magnetized rocks, such as oceanic sediments, very rapidly. Their 
principles of operation and various applications have been 
described by Goree and Fuller (1976). 

Although large quantities of data have now been obtained using 
this type of magnetometer, no detailed assessment of its operational 
characteristics has so far been published. Very little information 
has been provided by the manufacturers on such important details as 
alignment accuracy of the detector coils and the spatial uniformity 
of their response characteristics. In this paper I shall describe 
a series of tests run on an instrument in operation at NASA/Goddard 
Space Flight Center. This instrument is of similar design to almost 
all the superconducting magnetometers presently in use. 


I 


1 



DESCRIPTION OF THE INSTRUMENT 


Figure 1 Is a schematic showing the position of the single 
loop sense coils relative to the sample access, which in this 
instrument has a diameter of 6.8cm. The coils and sensors are 
surrounded by a superconducting lead shield which has the capability 
of "freezing in" the field in which it is allowed to cool below the 
shield transition temperature. 

It is clear from Figure 1 that only the vertical axis coil is 
of the true Helmholtz design, the two horizontal coils having a 
rectangular shape which in turn is curved around the dewar wall. 
Additional shielding is used above and below the detector coils to 
ensure the sample moment only couples with the coils within a limited 
region. In this way a zero reading can be taken with the sample 
only a small distance away from the detection region. The outputs 
of each coil sensor are calibrated by the manufacturers who use a 
small solenoid whose equivalent moment is known. 

Since each coil measures the total (remanent plus induced) moment, 
along each axis, it is necessary to make remanent moment measurements 
in a low field. How large a field can be tolerated depends on the 
magnetic characteristics of the samples being measured. In practice, 
a field smaller than 20 nT is usually found to be satisfactory. 

In this study, the maximum magnetic field within the detection region 
was less than 10 nT. The field gradients were less than 1 nT/cm 
(1 nT = 10” 5 gauss = 1 gamma). 
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RESPONSE OF THE DETECTOR COILS 


The spatial variation of the output of each of the three axes 
was checked using a small (1mm diameter) sample, which was given a 
saturation moment along one of its axes. Profiles of output versus 
vertical position were then obtained for each of the twelve horizontal 
positions shown in Figure 2. These profiles are shown in Figure 3- 
The sample was oriented so that its moment was along 
the axis being tested in each case. The figures adjacent to each 
profile represent the value of the output (in arbitrary units) at the 
midpoint of each profile. 

The results indicate that the region within which the response 
is uniform to within a few percent is fairly limited. The vertical 
axis defines the minimum volume of uniform response, which is a 
cylinder approximately 4 cm long and 2.5 cm in diameter, within which 
the variation in output is less than 5%. 

The horizontal axes responses are similarly uniform within a 
larger volume - a cylinder about 6 cm long and 3 cm in diameter. 
Although these volumes include the standard 2.5 cm diameter cores 
used in most paleomagnetic studies, it must be emphasized that these 
results are for a 6.8 cm access instrument. For the more common 
3.8 cm access instrument, the volume of uniform response would be 
correspondingly smaller, and almost certainly less than the standard 
core volume. The amount of error that this would cause in a single 
measurement would depend on the degree of inhomogeneity in the sample 
being measured. However, it is clear that significant errors could 
be introduced in this fashion. 
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CALIBRATION PROCEDURES : 

Absolute Calibrations were made using a small (0.476 cm diameter) 

5 turn coil. The equivalent dipole moment of this coil was taken as 

M = nIA/10 emu 
2 

where A is the area of the coil in cm , n is the number of turns and 
I ■is the ( current in amperes. 

The coil was positioned at the mid-point of the vertical response 
profile shown In Figure 3 (Z) and the current varied by an amount 
corresponding to 10"^ emu (1.12 mA). The calibration potentiometers 
in the magnetometer electronics were then adjusted to give corresponding 
changes in the magnetometer readouts of 10"^ emu, with the coil 
oriented along each axis in turn. The calibrations agreed to within 
a few percent with those performed by the manufacturers with the 
exception of the X sensor, which had recently been replaced and not sub- 
sequently recalibrated. 

Goree and Fuller (1976) have reported significant non-linearities 
in the output of their instrument. Linearity was therefore checked 
on the xlO range by plotting magnetometer output versus coil current, 
with the coil. The results indicated that linearity was better than 
0.5% for all three axes. The experiment was also repeated for the 
XI 00 range with similar results. 
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AXIS ORIENTATIONS 


The effective orientation of each axis was checked using a 
sample with most of its moment directed along a single axis. The 
direction of the major moment was then positioned at ten degree 
intervals in a plane at right angles to the axis being tested. In 
this way, cross-coupling effects show up as a variation in the 
output of the axis being tested, which should remain constant if 
the axes are truly orthogonal. The results are shown in Figure 4 
(a), (b), and (c). They indicate small amounts (about 1%) of 
crrss coupling between all three axes. Since the variations are 
approximately sinusoidal, it should be possible to correct for these 
effects using cross coupling coefficients whose values are known. 

If the true values of the three orthogonal components in the sample 
are X, Y, Z and the observed values x\ y\ l\ then 

X ’ ■ X + V + C zx Z 

yl - v + C xy X + C zy Z (2) 

z 1 • z + c xz x + c yz z (3) 

The cross-coupling terms C , C ... can be estimated from the 

Ajr At 

plots in Figure 4. However, a more convenient way to estimate them 

is to measure a calibration sample in different orientations. The 

calculations are considerably simplified if the calibration sample 

has most of its moment along one axis, as each cross-coupling term can 

then be estimated independently. Two measurements are then required 

for each cross coupling term. For example, to obtain C , the major 

yx 

moment of the sample is oriented along the Y axis, and the value 
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measured. The Y and Z components are then inverted by rotating the 

sample 180° about its X axis and a second value of obtained. 

It then follows |rom equation (1) that the difference between the 

first and second values of X 1 will be 2C Y, provided the last term, 

yx 

C 2X Z, is assumed negligible in both Instances. If Y 1 Is assumed 

equal to Y, the coefficient C can then be obtained. A more complete 

yx 

determination of the coefficients can be made using appropriate 
measurements on any sample at all. However, this will clearly 
Involve inversion of a matrix which the above method avoids. 

Shown in Table 1 are the results of measurements of a calibration 
sample with its cubical sample holder oriented in all possible (24) 
positions relative to the magnetometer measuring axes. This data 
can be used to obtain four independent estimates of each cross coupling 
term by suitably combining the data in the manner already described. 

The average values of the cross-coupling terms so obtained also 
appear in Table 1. The additional terms C xx , , and C zz are 
calibration normalisation terms. These were calculated by averaging 
all the measurements of the main component of magnetization for 
each axis, then normalising the results to the average value for 
the Z axis. In this way, calibration errors between the three axes 
are corrected for. The results of applying the corrections according 
to equations (1), (2), and (3) appear in Table 2. It is apparent 
that the scatter in the results has been reduced from almost 2° to less 
than 0.2°. Tables 3 and 4 show the results of similar measurements 
on a 2.5 cm diameter core of fine-grained basalt. Again, a 
considerable reduction in the scatter of the results is evident. 
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This reduction In scatter can be more easily seen In Figure 6, which 
shows plots of Inclination versus declination before and after correction 
for cross-coupling for the basalt sample. The slightly larger scatter 
of the corrected results In Table 4 (about 0.5°) Is probably due to 
the larger size of this sample compared to the original calibration 
sample, which was only ~ 2 mm in diameter and 2 mm In length. The 
directional variation may therefore be due to small Inhomogeneities 
in the rock sample causing differences in the coil responses as 
outlined in the previous section. 

One interesting feature of the cross coupling coefficients shown 

in Table 1 is that, for example, C is not the same as C . This 

xy yx 

implies that the coupling cannot simply be attributed to misalignment 
of the coils, for this would clearly imply that these two terms should 
be equal. Some of the coupling is probably between the wires leading 
up to each coil and che other coils. This would explain the no~- 
commutative behavior of the coupling terms. 


NOISE MEASUREMENTS 


Shown In Figure 5 are output traces recorded over a 5 minute 

* 

period using the Ihz bandwidth filter. The traces also show deflections 

recorded for Insertion of samples having the indicated moments. The 

-8 

short term noise varies from about 2x10 emu r.m.s. for the x axis 
-8 

to 1x10 emu r.m.s. for the Z axis. The long term drift ranoes from 

•9 -8 

5x10 emu/min for the Z axis to 2x10 emu/min for the V axis. The 

larger drift in the horizontal sensor outputs seemed to be related to 

less effective shielding in the horizontal directions. It was found 

experimentally that a magnet producing 10 gauss at the sense coils 

-5 

caused a change of 10 emu in the outputs of the horizontal sensors. 
This implied that field changes of 0.01 gauss, which are fairly common 
in urban evnironments#can cause observable noise in these magneto- 
meters outputs. The vertical shielding factor is much greater - the 
10 gauss magnet produced no observable change in the output of the 
vertical sensor. 


DISCUSSION 


Superconducting magnetometers of this type are clearly useful 
instruments for measuring weak magnetization in paleomagnetic samples. 
Although their spatial response characteristics are clearly more 
uniform than those of either spinner magnetometer coils or fluxgates, 
significant variations are nevertheless present. 

To eliminate the effect of these, measurements need to be 
performed with samples in different orientations and the results 
averaged as has been done with spinner magnetometer data. It has 
been shown that with suitable precautions, a ore- inch diameter core 
can be measured with a single insertion to an angular accuracy of 
better than 0.5°. However, with a 3.8 cm access Instrument, the angular 
error could be much greater. 
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TABLES 


Table 1 


Table 2 


Table 3 


Table 4 


- Measurements of calibration sample in different orientations. 
Also shown are the cross coupling terms calculated from this 
set of measurements using the methods described in the text. 

- Corrected calibration sample measurements. Note the large 
reduction in scatter compared to Table 1. 

- Measurements similar to those in Table 1, for a 2.5 cm dia- 
meter core of fine grained basalt. 

- Corrected results for the basalt core. 



FIGURE CAPTIONS 


Figure 1 - Sense coil geometry for a two-axis system. 

Figure 2 - Horizontal positions of response profiles. The numbers in 
the circles correspond to the numbers on the left hand side 
of each of the profiles in Figure 3. 

Figure 3 - Response of each axis to a sample as a function of vertical 
position in each of the positions defined in Figure 2. The 
numbers at the center of each profile are the output moment 
values (arbitrary units) at the optimum position defined by 
the vertical responses. (See text) 

Figure 4 - Cross coupling response for: (a) the X axis, (b) the Y axis 

and (c) the Z axis. The upper two curves in each case are the 

variations in output of two of the axes when the sample is 

rotated about the third axis. The lower curve is the varia- 
tion in the response of this third axis, which for no cross 

coupling, should be zero. 

Figure 5 - Noise measurements made over 5 to 7 minute periods (not 

synchronously) for each axis. Also shown is the effect of 
inserting the small moments indicated, into the sense region. 

Figure 6 - Plots of inclination, versus declination for the basalt core 

as measured in all 24 orientations, before and after correction 
for cross-coupling effects. 
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Table 2: Corrected measurements of Calibration Sample 
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Table 3: Uncorrected measurements on 2-5 cm core of fine grained basalt 
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Table 4: Corrected results for same sample as Table 


Sensor Output 



Fig. I SENSE COIL GEOMETRY (TWO-AXIS SYSTEM) 
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